A Direct Probe of Electronic Nematic Order; Symmetry Information in Scanning 

Tunneling Microscope Images 
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An electronic nematic state spontaneously breaks a point-group symmetry of an underlying lattice. 
As a result, the nematic- isotropic transition accompanies a Fermi surface distortion. However, the 
anisotropic nature of the nematic state at a macroscopic scale can be easily wiped out when domains 
of different orientations of nematic order exist. We suggest that a spatial pattern of local density 
of states (LDOS) in the presence of a non-magnetic impurity can be a direct probe of the nematic 
order. We study various patterns of LDOS across the quantum phase transition between the isotropic 
and nematic phases. Especially the Fourier transformed local density of states (FT-LDOS), which 
can be deduced from scanning tunneling microscope images, represent a transparent symmetry of 
an electronic structure. The application of our results to the bilayer ruthenate, Sr3Ru207 is also 
discussed. 

PACS numbers: 71.10.Hf,71.20.-b,71.55.-i,73.22.Gk,73.43.Nq 



Introduction: Recently, the existence of an intermedi- 
ate form of matter dubbed the electronic nematic liquid 
was proposed to explain novel features observed in ultra- 
pure Sr3Ru207.[H Especially, a large magnetoresistive 
anisotropy in the close vicinity of a metamagnetic quan- 
tum critical point shows a striking similarity to those 
reported in ultra-clean two-dimensional electron systems 
in high magnetic fields [2|, u3] , where the electronic nematic 
phase was also suggested. The metamagnetic transition 
in SrsRu207 possesses a critical end point at a finite tem- 
perature, and it was proposed that a quantum critical 
point can be achieved by tilting the magnetic field to- 
wards the direction perpendicular to the Ru-0 plane. Q 
However, when the system gets close to the putative 
quantum critical point, a new phase bounded by two suc- 
cessive metamagnetic transitions was found in an ultra- 
clean sample. [5] The emergent phase is characterized by 
a high residual resistivity compared to the phases nearby, 
and it was suggested that a formation of the electronic 
nematic phase leads to these remarkable phenomena of 
consecutive metamagnetic transitions and the high resid- 
ual resistivity. 0, @, 0] 

The electronic nematic phase is characterized as an 
anisotropic metallic state with a broken point-group sym- 
metry of an underlying lattice. SlESIHEIi! For ex- 
ample, on a square lattice, there are two distinct nematic 
phases depending on how the 7r/2 rotational symmetry 
is broken. [14] A consequence of the broken symmetry is 
a spontaneous Fermi surface distortion, which leads to 
an anisotropy between two directions of the longitudinal 
resistivity. Indeed, a large magnetoresistive anisotropy 
was observed in Sr3Ru207 consistent with a formation 
of nematic phase. On the other hand, when the field is 
perpendicular to the Ru-0 plane, the anisotropy is no 
longer present, but the high residual resistivity becomes 
pronounced. It was proposed that the nematic domain 
formation results in the high resistivity in the emergent 



phase, which involves two different orientations of the ne- 
matic order in the square lattice system. [7[ In the pres- 
ence of the nematic domains, the nature of the broken ro- 
tational symmetry can be hardly observable in any bulk 
measurement. 

In this paper, we suggest that a spatial pattern of the 
LDOS as measured by scanning tunneling microscopy 
(STM) provides a direct probe of the nematic order. Note 
that the nematic-isotropic transition is first order, thus 
the insignificant effect of order parameter fluctuation on 
the Fermi surface measurement is in contrast to the case 
near the nematic quantum critical point. [15[ We show 
that the Friedel oscillations around an impurity show 
visible anisotropic patterns between x and y directions 
in the nematic phase, as one expected. In particular, the 
FT-LDOS shows a transparent difference between the ne- 
matic and isotropic phases across the transition by the 
magnetic field. Therefore, STM on the bi-layer ruthen- 
ate is highly desirable for the confirmation of the nematic 
phase emerging near the vicinity of the putative quantum 
critical point. 

Phenomenological Model: The essential physics of the 
broken symmetry possessed in LDOS is insensitive to the 
choice of a Hamiltonian, as long as the nematic-isotropic 
transition is captured in the Hamiltonian. On the other 
hand, the detailed patterns of the LDOS depend on the 
interactions as well as band structure. To capture an elec- 
tronic structure symmetry embedded in spatial patterns 
of LDOS across the transition, we use a phenomenolog- 
ical Hamiltonian with a quadrupole interaction flpj , 
which successfully described the phenomena of the meta- 
magnetic transitions, high resistivity, and anisotropic 
magnetoresistivity in the bilayer ruthenate. 0, 0, 

Since the nematic phase contains a broken rotational 
symmetry, an analog of the order parameter of the ne- 
matic liquid crystal is the first step towards the construc- 
tion of a phenomenological model. H, 12, 16[ We defined 
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FIG. 1: Mean field result of nematic order (solid line) and 
magnetization (dashed line) as a function of magnetic field. 
The points marked as (a), (b), (c), and (d) are where we 
present LDOS in the following figures. The inset shows the 
Fermi surface of up-spin and two dominant extreme vectors 
in the nematic phase discussed in the main text. 



Q(q), a quadrupole tensor on the square lattice as 



,(q) = £' 
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where ("l (k) = cos k x — cos k y and ("2 (k) = 2 sin k x sin k y . 
We can write down the Hamiltonian as follows using 
the quadrupole tensor, which leads to the nematic order 
within the mean-field theory. 

H =^( ek - M ^^) c t (j c k(J -^F 2 (q)TV Qt (q)4,(q), 



qcr 



(2) 

where ek = —t(cosk x + cosk y ) is the electronic dis- 
persion of the tight-binding model on the square lat- 
tice, and ji is a chemical potential. o\±bB is a Zee- 
man term, where cr = +1 for up-spin and a = —1 
for down-spin, ^(q) represents the strength of the 
quadrupole-quadrupole interaction, which can be viewed 
as the Landu quasiparticle interaction of I = 2 angular 
momentum channel. 0, [H, [3, 0, EH 

There are two types of the Fermi surface distortion as- 
sociated with £1 and (2 which represent the Fermi surface 
elongations along the directions parallel to the axis and 
diagonal to the axis, respectively. The previous mean- 
field results showed that for a given value of F2, the 
diagonal distortion does not occur, so we focus on the 
nematic phase with a Fermi surface distorted along the x 
or axis of the square lattice. The mean field equation 
for the uniform nematic order can be written as 
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where the nematic order parameter is defined as A a = 
^2(0X1 (k) ^c kcr c kcr ^ . The order parameter as a func- 
tion of magnetic field H for a given chemical potential 



was studied in RefH The strongly first order nature 
of transition is originated from the van Hove singular- 
ity, and the density of up (or down)-spin as a function 
of Zeeman field jumps at the transitions, which results 
in two consecutive metamagnetic transitions. In Fig. HJ 
we showed the earlier result of the order parameter to 
indicate the locations in the phase diagram where we 
present LDOS. We set F 2 N = 0.1 and /iN = -0.02, 
where No = 2ir 2 t. (a) represents the isotropic phase 
with ATVo = 0, and (b)-(d) the nematic phase with 
A7V = 0.0113,0.0116 and 0.0126, respectively. We car- 
ried out the standard T-matrix formalism to compute 
FT-LDOS and show the results for (a) -(d) cases below. 

Local density of states with a single impurity: The 
FT-LDOS, n(q;cj) with a single impurity is obtained us- 
ing the standard T-matrix formalism. 

n(q;^) = n (q;^)+Im J dkG^ (k+q; w)T aa ' (w)G$ (k; uj), 

(4) 

with the T-matrix, T _1 (cj) = V' 1 - Jdk G^ 0) (k;cj), 
where V s represents the strength of the impurity, and 
we set V s No = —0.005. The single particle retarded 
Green's function, modified by the nematic order self- 
consistently can be written as follows. 
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Here we ignore effects of the impurity on the nematic or- 
der, so our results are valid only when the strength of an 
impurity and the bias voltage is smaller than the energy 
scale of nematic order. We use 2-dimensional adaptive 
integration to compute n(q; w) shown in Fig. [2j 

Fig. [2] shows LDOS and FT-LDOS at the Fermi en- 
ergy {u = 0) for various strengths of nematic order cor- 
responding to (a)-(d) indicated in Fig. [TJ The left three 
columns show LDOS in real space for up-spin, down-spin, 
and sum of both spins from left to right, respectively. FT- 
LDOS images are also shown in the next three columns. 
Each row corresponds to (a) to (d) cases in Fig. [T] from 
top to bottom. The color represents a relative strength of 
LDOS (FT-LDOS) induced by the impurity. The Friedel 
oscillations around the impurity located at the center 
show visible anisotropic patterns between x and y direc- 
tions. The peak structures in FT-LDOS are associated 
with the Fermi surface nesting vectors. In the isotropic 
phase of (a), the peaks of the FT-LDOS form an aster- 
oidal shape with the strongest peak near (±7r,±7r) which 
is the nesting vector of the half-filled tight-binding dis- 
persion. In the nematic states, (b) and (c), the asteroidal 
peak structure abruptly changes to a sandglass-like shape 
due to the opening of the Fermi surface along the x-axis. 
The anisotropic shape in the nematic phase is mainly due 
to the extreme vectors of 2kpi and 2k^2 as shown in the 
inset of Fig. [lj These extreme vectors also determine 
the wave length of the Friedel oscillations which will be 
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FIG. 2: LDOS and FT-LDOS at oo = for the marked points, (a) - (d) in Fig. [Q The left three columns are LDOS and the 
right tree columns are FT-LDOS. Each set of three columns shows spin- up, spin-down, and sum of both spins from left to right. 
From top to bottom, nematic order increases as marked in Fig. [1] 



further discussed below. While the nematic order is max- 
imized in (d), we found that both FT-LDOS and LDOS 
do not show a significant signal of the nematic phase be- 
cause of the destructive interference between two extreme 
vectors of 2kp2 and 2kpi. However, the nematic order 
signals at finite bias voltages represented by the dashed 
and dotted lines in the bottom of Fig. [3l 

FT-LDOS along q x - and g y -axis for finite uo are shown 
in Fig. [3j The dramatic changes in FT-LDOS occur 
across the transition. FT-LDOS along the q x direc- 
tion lose their weight, while additional structures de- 
velop along the q y direction. One is originated from 
the extreme vector of 2kp2 connecting (7r, dzfe^), while 
the other from the extreme vector of 2kpi connecting 
(0, ±fcpi), of the Fermi surface shown in Fig. dJ These 
peaks along the q y direction give a unique pattern of 
Friedel oscillation in LDOS. LDOS along x and ?/-axis 
shown in Fig. 0] also reflect the broken symmetry in the 
nematic state. For example, LDOS represented by the 
dashed line clearly show different amplitude variations 
along the ^-direction, which roughly correspond to the 
wave length (kpi — kp2)~ 1 • 

The magnetic field dependence of LDOS is shown in 
Fig. [51 At the isotropic-nematic transition occurring 
around [ibH = 0.13t, LDOS at impurity site shows a 
sudden decrease (increase) for an attractive (repulsive) 
impurity potential (not shown here). It is interesting note 



that LDOS at neighboring sites such as (0, 1) (in unit of 
the lattice spacing a = 1), show a dramatic change across 
the transition, while LDOS at (1,0) does not show any 
significant change across the transition. The field de- 
pendence of LDOS also possesses the strong anisotropic 
nature of the nematic phase. 

Discussion and Summary: Strongly correlated elec- 
tron systems with d or /-electrons have demonstrated 
rich phase diagrams with many different phases. When 
the transition temperature between two different phases 
is driven to zero by the application of a magnetic field, 
pressure, or via chemical doping, we can achieve the 
quantum critical point (QCP). The behavior of the sys- 
tem close to the QCP is very different from the behavior 
in the vicinity of the classical phase transition, and it 
has been an intensive research area. Along the journey 
of searching for QCP, we often discovered that a novel 
phase emerges at the vicinity of the desired QCP. 

The emergent phase at the vicinity of the QCP in the 
bilayer ruthenate is characterized by the anisotropic mag- 
netoresistivity, when the field is tilted from c-axis. On the 
other hand, the high resistivity without anisotropy was 
reported, when the field is along c-axis, H || c. It was pro- 
posed that the nematic domain formation is responsible 
for the high resistivity for H || c. These imply that the 
nematic domain may align when the field is tilted from 
c-axis mostly due to the in-plane component of magnetic 
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FIG. 3: FT-LDOS (subtracted by FT-LDOS in the absence of 
an impurity) along q x - and q y -axis are shown for the marked 
points (a), (b), and (d) in Fig. Q] The left-half is FT-LDOS 
along g^-axis and the right- half, along q x -axis. The red solid 
line, the magenta dotted, and the green dashed lines corre- 
spond to FT-LDOS at co = 0, —O.lt, and O.lt, respectively. 
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FIG. 5: LDOS (at uj = 0) vs. magnetic field at several neigh- 
boring sites (with lattice spacing a = 1) around the impurity. 



nonmagnetic impurity across the nematic-isotropic tran- 
sition tuned by the external magnetic field, and found 
the followings. (1) FT-LDOS as well as the Friedel oscil- 
lations in real space show a transparent indication of the 
broken rotational symmetry. In particular, the Fourier 
transformed image of LDOS which can be deduced from 
STM data gives direct information on the topology of a 
Fermi surface associated with the electronic nematic or- 
der. (2) Across the transition, the change of the Fermi 
surface volume accompanies a sudden change in LDOS. 
Therefore, LDOS at the impurity and neighboring sites 
show significant and anisotropic changes across the tran- 
sition. (3) In the nematic phase, it is possible that the 
anisotropic signal can be washed out, when the signals 
from two extreme vectors have destructive interference. 
However, anisotropic properties are still detectable at fi- 
nite UJ. 

We thank J. C. Davis for encouraging us to present 
the current work. This work was supported by NSERC 
of Canada, Canada Research Chair, Canadian Institute 
for Advanced Research, and Alfred P. Sloan Foundation. 
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FIG. 4: LDOS for u — are shown along x- and y- axis with 
an impurity at the origin. The left half of the figure shows 
LDOS along y-axis, and the right half, along x-axis. The 
solid and dashed lines correspond to ANo = and 0.0116, 
respectively. 

field, and results in the anisotropic magnetoresistivity. 
At present, we do not have a microscopic picture which 
describes interplay between the nematic domain forma- 
tion and the direction of magnetic field, which needs to 
be further investigated in future. 

Here we proposed that a local probe such as STM is in- 
sensitive to a macroscopic inhomogeneity, and therefore 
it can directly confirm an existence of the nematic order. 
We studied FT-LDOS and LDOS in the presence of a 
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